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A model of coherent tunneling, which combines multi-orbital tight-binding approximation with 
Landauer-Biittiker formalism, is developed and applied to all-semiconductor heterostructures con- 
taining (Ga,Mn)As ferromagnetic layers. A comparison of theoretical predictions and experimental 
results on spin-dependent Zener tunneling, tunneling magnetoresistance (TMR), and anisotropic 
magnetoresistance (TAMR) is presented. The dependence of spin current on carrier density, mag- 
netization orientation, strain, voltage bias, and spacer thickness is examined theoretically in order 
to optimize device design and performance. 

PACS numbers: 75.50.Pp, 72.25.Hg, 73.40.Gk 



I. INTRODUCTION 

Metallic magnetic tunnel junctions (MTJs) are impor- 
tant building blocks of the already existing spintronic 
devices, such as magnetic random access memories, mag- 
netic heads, and sensors. The MTJ structure consists of 
ferromagnetic layers separated by a thin insulating bar- 
rier through which the carriers tunnel. The resistance of 
such a junction depends on the relative alignment of the 
magnetization vectors in the ferromagnetic layers, i. e., 
the structure exhibits the tunnelling magnetoresistance 
(TMR) effect. 

All-semiconductor MTJ structures offer potential for 
precise control of interfaces and barrier properties, par- 
ticularly in the case of III-V compounds, for which epi- 
taxial growth of complex heterostructures containing fer- 
romagnetic (Ga,Mn)As or (In,Mn)As layers is especialy 
advancedi Extensive studies of MTJs with (Ga,Mn)As 
ferromagnetic contacts carried out by various groups re- 
sulted in an increase of the observed TMR ratio from 
about 70% reported by Tanaka and Higc 2 to values 
higher than 250%A^ Another key factor for develop- 
ing novel functional semiconductor spintronic devices is 
an efficient electrical injection of spin polarized carriers. 
Here again the III-V ferromagnetic p-type semiconductor 
(Ga,Mn)As with its high spin polarization^, appears as a 
promising material. The electrical spin injection from 
p-(Ga,Mn)As into non-magnetic semiconductor has first 
been achieved by injection of spin polarized holes J Later, 
injection of spin polarized electrons was demonstrated 
employing interband tunneling from the valence band of 
(Ga,Mn)As into the conduction band of an adjacent n- 
GaAs in a Zener-Esaki diode£ Recently, a very high spin 



polarization of the injected electron current (ca 80%) was 
obtained in such devicesAifi It should be mentioned that 
both effects, the TMR and the spin polarization of tun- 
neling current in the Zener-Esaki diode, decrease rapidly 
with the increase of the applied bias - a phenomenon ob- 
served also in the metallic TMR structures, and still not 
fully understood. Finally, it seems that TMR is sensitive 
to the direction of the applied magnetic field in respect to 
the direction of current and crystallographic axes. This 
so-called tunnel anisotropic magnetoresistance (TAMR) 
effect was observed in structures containing a single fer- 
romagnetic electrod o 11 ! 12 as well as in typical TMR MTJ 
with two ferromagnetic contacts^i^ 

These challenging experimental findings call for a the- 
ory that would describe the tunneling in semiconductor 
MTJs and would indicate the ways for optimized de- 
sign of the devices. Since the ferromagnetic coupling 
in (Ga,Mn)As is mediated by the holesji*& a meaning- 
ful theory has to take into account the entire complex- 
ity of the valence band, including the spin-orbit inter- 
action. Furthermore, the intermixing of valence bands 
caused by spin-orbit coupling shortens the spin diffu- 
sion length and makes it comparable to the phase co- 
herence length. This renders the models based on the 
classical spin-diffusion equation, which describe satisfac- 
torily spin transport phenomena in metallic MTJs, non 
applicable directly to the structures containing layers of 
hole-controlled diluted ferromagnetic semiconductors. 

The model of vertical transport in modulated struc- 
tures of magnetic semiconductors described here com- 
bines the two-terminal Landauer-Biittiker formalism 
with the empirical multi-orbital tight-binding descrip- 
tion of the semiconductor band structure. In this way, 
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the quantum character of spin transport over the length 
scale relevant for the devices in question is taken into 
account. Furthermore, the tight-binding approach, in 
contrast to kp models employed so-far^ allows for a 
proper description of effects crucial for spin transport in 
heterostructures such as atomic structure of interfaces, 
effects of Rashba and Dresselhaus terms as well as tun- 
neling involving k states away from the center of the 
Brillouin zone. Our model has recently been applied to 
describe selected features of Zener-Esaki diodesiSiii and 
TMR devices^ as well as it was adopted to examine an 
intrinsic domain-wall resistance in (Ga,Mn)Asii& 

The remaining part of the present paper is organized 
as follows. In Section II we present our model, specifying 
the tight-binding parametrization, scattering formalism, 
and the way transfer coefficients are determined. In Sec- 
tion III, the calculated dependencies of spin polarization 
of the current in the Zener-Esaki diode on carrier den- 
sity, trigonal distortion, and magnetization direction are 
shown and compared with available experimental find- 
ings. The calculated TMR and TAMR ratios for struc- 
tures containing two ferromagnetic (Ga,Mn)As contacts, 
and their dependencies on the voltage bias, trigonal de- 
formation, and on the width of non-magnetic spacer layer 
are presented and discussed in Section IV in reference 
to experimental results. Section V contains conclusions 
emerging from our work, particularly concerning possible 
optimizations of device performance. 



II. THEORETICAL MODEL 

We consider a prototype heterostructure, which is uni- 
form and infinite in the x and y directions and has mod- 
ulated magnetization along the z growth direction. The 
heterostructure is connected to two semi-infinite bulk 
contacts denoted by L and R and biased. In all cases 
considered, spin polarized carriers are injected from the 
ferromagnetic left lead. Our goal is to calculate the elec- 
tric current in the structure and the degree of current 
spin polarization outside the left lead. Typical length 
of the studied structures is comparable to the phase co- 
herence length. Therefore, we restrict ourselves to the 
vertical coherent transport regime that we treat within 
Landauer-Buttiker formalism, where the current is deter- 
mined by the transmission probability from the ingoing 
Bloch state at the left contact to the outgoing Bloch state 
at the right contact. In the presence of spin-orbit cou- 
pling, the spin is not a good quantum number. The only 
preserved quantities in tunneling are the energy E and, 
due to spatial in-plane symmetry of our structures, the 
in-plane wave vector feii. We use semi-empirical tight- 
binding formalism to calculate electronic states of the 
system for given fc|| and E and further to compute trans- 
mission coefficients. 



A. Tight-binding model 

First, we describe the construction of the tight-binding 
Hamiltonian matrix for 'normal' GaAs and AlAs as well 
as ferromagnetic (Ga,Mn)As layers of the heterostruc- 
ture. To describe the band structure of the bulk GaAs 
and bulk AlAs, we use the nearest neighbor (NN) sp 3 d 5 s* 
tight-binding Hamiltonian (resulting in 20 spin-orbitals 
for each anion or cation), with the spin-orbit coupling 
included. 19 This model reproduces correctly the effective 
masses and the band structure of GaAs and AlAs in the 
whole Brillouin zone. With the tight-binding Hamilto- 
nian introduced above, each double layer (cation + an- 
ion) is represented by 40 x 40 matrix. It should be pointed 
out that the d orbitals used in our sp 3 d 5 s* parametriza- 
tion are not related to the 3d semi-core states and are 
of no use for description of Mn ions incorporated into 
GaAs. The presence of Mn ions in (Ga,Mn)As is taken 
into account by including the sp-d exchange interactions 
within the virtual-crystal and mean-field approximations. 
In the spirit of the tight-binding method, the effects of 
an external interaction are included in the on-site di- 
agonal matrix elements of the tight-binding Hamilto- 
nian. Here, the shifts of on-site energies caused by the 
sp-d exchange interaction are parameterized in such a 
way that they reproduce experimentally obtained spin 
splitting: Noa — 0.2 eV of the conduction band and 
N /3 = -1.2 eV of the valence band^S It should be, 
however, mentioned that since we neglect exchange in- 
teractions between the holes, the valence spin splitting 
is presumably underestimated by about 20% & The other 
parameters of the model for the (Ga,Mn)As material and 
for the NN interactions between GaAs and (Ga,Mn)As 
are taken to be the same as for GaAs. This is well moti- 
vated because the valence-band structure of (Ga,Mn)As 
with small fraction of Mn has been shown to be quite sim- 
ilar to that of GaAs«2& Consequently, the valence band 
offset between (Ga,Mn)As and GaAs originates only from 
the spin splitting of the bands in (Ga,Mn)As. The Fermi 
energy in the constituent materials is determined by the 
assumed carrier concentration and is calculated from the 
density of states obtained for tight-binding Hamiltonian. 
Our calculations of the Fermi energy for various hole con- 
centrations are consistent with the corresponding results 
presented in Ref . |fl Having determined the Hamiltonian 
of the system, we are now in the position to define the 
current and current spin polarization in the presence of 
spin-orbit coupling. 



B. Current and current spin polarization 

For a given energy E and in-plane wave- vector /c|| , the 
Bloch states in the left L and right R leads (i and j, 
respectively) are characterized by the wave vector com- 
ponent k± perpendicular to the layers and are denoted 
by |i, and \R, kji ±j}, respectively. The indices i 

and j indicate all possible pairs for 40 bands described by 
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the tight-binding Hamiltonian. The transmission proba- 
bility Tl jz l ± t—tR^R ±j is a function of the transmission 
amplitude tL,k L _!_,»— >--K,feB ± a (E, fell) and group velocities 
in the left and right lead, fz,,_L,i and VR t ±j 

I 2 «fl,-L,j 



Fi,fcl,,i, ! ^-R,fcfl,_L, 3 ( £, > fc || 



The current flowing in the right direction can now be 
written as2i 



d 2 k«dEf L (E) 



{2-Kfh JBZ 



(2) 



VL,±,i,v R ,±.j>0 

where Jl or respectively Jr are the electron Fermi dis- 
tributions in the left and right interface and i, j number 
the corresponding Bloch states. Plugging in the expre- 
sion given in Eq. 1 and using the time reversal symmetry 



T L ,k L:± , i ^R,k R , ± jE,k ll )= (3) 

= T L-k Rt±tJ ^R,-k L , ±tZ (E, fe|| ) 

we get 

— e 



J 



d 2 k«dE[f L (E)-f R (E)] 



(4) 



BZ 



iiL,fc L ,_ L , i ^_R,fe Ji ,_ LJ (^, fen) 



VL,±,i,VR, ±,j>0 



|2 VR,_L,j 
VL,±.i ' 



Let us define the spin polarization of the outgoing 
Bloch state, in respect to magnetization direction in the 
source lead 

P R ,kR^{E,k\\) = (R^R.^in-slR^kR.^), (5) 

where fi is the magnetization direction vector and s is 
the spin operator. Then, we can define the spin polarized 
current 
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:— 6 — f d 2 k^dE[f L {E)-fR{E)] 

[2Tryh J BZ 

^2 T L.k Li±il ^R : k R _ ±J (E,k l \)PR tkR ± 



(6) 



VL,±,i,VR,±.j>0 



The spin polarization of the coherently transmitted cur- 
rent is now equal to 



P, = 



3 



(7) 



To calculate the current one has to determine the 
transmission probability, thus the transmission ampli- 
tude ti, t k L J _ li ->fe R R x j (E, fen) and the group velocities 



v L,±.j of the ingoing and Ur,_lj of outgoing states. These 
can be obtained by solving the Schrodinger equation for 
the structure with the appropriate scattering boundary 
conditions. In our studies, we follow closely the proce- 
dure detailed in Refs. [52 and I23L which we have gener- 
alized to the case with spin-orbit coupling. 



III. SPIN-DEPENDENT ZENER TUNNELING 



Effect of carrier densities 



The rather high « 80% spin polarization of the tunnel- 
ing current, obtained recently in Zener-Esaki diodesSiifi 
opens new perspectives for applications of electron spin 
injection. The degree of current spin polarization de- 
creases sharply with the biasj2ii£ an effect explained 
quantitatively by our models On the other hand, it is 
well known that magnetic characteristics of (Ga,Mn)As 
depend strongly on both hole and manganese concen- 
trations and that (Ga,Mn)As films exhibit a variety of 
anisotropic properties It is thus obvious that the 
degree of spin polarization of the tunneling current may 
depend on these intrinsic features of (Ga,Mn)As layers. 
Indeed, we have already showni£ that a higher content 
of magnetic ions x in Gai-^Mn^As results in an increase 
of the spin polarization of the tunneling current. In con- 
trast, an opposite change was obtained when the hole 
concentration was increased. 

In order to get a better insight into processes control- 
ling spin polarization of the current, we have examined 
the dependence of tunneling on the in-plane wave vector 
in the low bias limit. As shown in Fig. ^a), the total 
current is dominated by the tunneling from states close 
to the r point. This is because in the tunneling pro- 
cess the in-plane components of fc|| wave vectors are con- 
served, i. e., they have to match to the small fc|| vectors 
at the Fermi level in the conduction band of n-type GaAs 
(n = 10 19 cm -3 ). Turning to current spin polarization, 
we note that it decreases with the hole concentration p 
because the higher p the smaller is the spin polarization 
at the Fermi level in the vicinity of the center of the Bril- 
louin zone. This is shown in Fig. [2] where the cross sec- 
tions of the Fermi sphere for different hole concentrations 
in p-Gao.92Mno.08As with the saturated value of magne- 
tization are presented. At the same time we find that 
the large fc|| vectors are responsible for the spin polariza- 
tion of the current, as depicted in Fig. db). This sug- 
gests that a higher concentration of electrons in n-GaAs 
layer should result in matching of the larger k vectors 
and thus in higher spin polarization of the current. In- 
deed, when the electron concentration is increased from 
n = 10 19 cm" 3 to n = 10 20 cm" 3 , the current spin po- 



larization becomes higher by about 



4 



a) 




b) 



-0.09 
1-0.18 
- H - 0.27 - 
-0.36 
-0.45 



-0.73 




--0.4 



-0.4 -0.2 0.2 0.4 
k 



-0.4 -0.2 0.2 0.4 
k 



FIG. 1: [color on-line] Dependence of tunneling current (a) 
and its spin polarization (b) on the in-plane wave vector 
in p-Gai-ajMn^ As/n-GaAs Zener-Esaki diode in the limit 
of low bias. The calculation was performed for the hole 
concentration p = 3.5 x 10 20 cm" 3 , electron concentration 
n = 10 19 cm" 3 , and the saturated magnetization correspond- 
ing to the Mn content x = 0.08. 
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FIG. 2: [color on-line] Cross section of the valence bands 
at the Fermi energy for various hole concentrations p in p- 
Gao.92Mno.08As. Color scale denotes spin polarization. Mag- 
netization is taken along the [110] direction. 



B. Anisotropic Zener tunneling - in-plane 
magnetization 

In this and next subsection we asses the importance 
of a new mechanism called tunneling anisotropic mag- 
netoresistance (TAMR). This effect consists of a change 
in the tunnel resistance upon the rotation of magneti- 
zation. The phenomenon, recently discovered in struc- 
tures with a single (Ga,Mn)As ferromagnetic layer^iiiS 
results from the fact that tunneling resistance depends 
on the relative orientation of magnetization in respect to 
the direction of current and crystallographic axes. This 
is due to the strong spin-orbit coupling and the highly 
anisotropic Fermi surface in (Ga,Mn)As (compare Fig. EJ) 
- it is why TAMR was not reported for structures based 
on ferromagnetic metals, where typically spin-orbit char- 



acteristic energies are smaller than the Fermi energy. 

We consider a simple junction consisting of several 
layers of p-type Gai_ x Mn x As, x = 0.08, p = 3.5 x 
10 20 cm -3 , followed by several layers of n-type GaAs, 
n = 10 19 cm" 3 in the weak bias limit. Interestingly, as 
shown in Fig. [3] the model reveals that the current mag- 
nitude and its spin polarization differ for magnetization 
along [110] and along [TlO] crystallographic axis, even 
in the absence of any extrinsic deformation. This re- 
flects the asymmetry of the (Ga,Mn)As/GaAs interface, 
at which the [110] and [TlO] directions are not equiva- 
lent, so that the symmetry of the zinc-blende crystal 
is reduced to Ci v for the heterostructure in question. Ac- 
tually, a 6% difference in spin polarization of the current 
for the unstrained structure, which is visible in Fig. EJb) , 
agrees with that observed experimentally^ 

The intrinsic anisotropy of the Zener tunneling cur- 
rent depends on the hole concentration, as shown in 
Fig. 01 The change in the tunneling current upon ro- 
tation of the magnetization vector from [110] to [TlO] 
direction increases with the decrease of hole concentra- 
tion, reaching 8.5% for p = 1 x 10 19 cm" 3 . Such a low 
hole concentration can, in fact, correspond to a deple- 
tion region in the p-(Ga,Mn)As/n-GaAs junctionjiii^ 
though hole localization may render our theory invalid 
in this low hole concentration range. Thus, our model 
predicts for (Ga,Mn)As/GaAs Zener diodes an in-plane 
TAMR magnitude of the order of several percents with- 
out assuming any extrinsic strain. Although the ob- 
tained TAMR value agrees with the observation reported 
in Ref. 111! it should be emphasized that in that ex- 
periment (Ga,Mn)As/A10 x /Au tunnel junction was ex- 
amined. Moreover, the symmetry of the experimental 
TAMR effect implies the existence of an extrinsic defor- 
mation breaking the equivalence of [100] and [010] crys- 
tallographic axes. 
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FIG. 3: [color on-line] Dependence of Zener tunneling current 
(a) and its spin polarization (b) on the direction of in-plane 
magnetization without strain and for deformations up to 1% 
applied along the [110] axis. Calculations performed for Mn 
content x = 0.08, hole concentration p — 3.5 x 10 20 cm -3 , 
and electron concentration n = 10 19 cm -3 . 



Typically (Ga,Mn)As films exhibit uniaxial anisotropy, 
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whose character implies the presence of an extrinsic trig- 
onal distortion along the [110] axisi2£ A strain as small as 
0.05% was found to explain the magnitude of the corre- 
sponding uniaxial in-plane anisotropy fields The effect 
of the trigonal strain on the Zener current and its spin po- 
larization is presented in Fig. together with previously 
discussed results for unstrained structures. As seen, the 
strain causes an additional in-plane anisotropy. How- 
ever, a rather strong deformation is needed to obtain a 
significant dependence of spin current polarization on the 
direction of the magnetization vector. Even for 0.1% de- 
formation, which is two times larger than that evaluated 
in Ref. [2|j, the anisotropy of the spin polarization of the 
current is still very small. When a strong, 1%, deforma- 
tion is assumed, the calculation predicts a 10% increase 
of current spin polarization for magnetization along [110] 
axis. At the same time, the obtained spin polarization 
for the magnetization along [110] direction is smaller by 
about 30%. 
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FIG. 4: [color on-line] Dependence of Zener tunneling current 
(a) and its spin polarization (b) on the direction of in-plane 
magnetization for various hole concentrations; Mn content 
x = 0.08 and electron concentration n = 10 19 cm -3 . 



C. Anisotropic Zener tunneling - out-of-plane 
magnetization 

The spin dependent interband tunneling is sensitive 
not only to the in- plan e magnetization direction. It has 
been shown in Ref. [12J that rotation of magnetization by 
applying an out-of-plane magnetic field leads also to a 
TAMR signal in a Zcncr-Esaki diode. The magnitude of 
perpendicular tunneling anisotropic magnetoresistance is 
defined as 



TAMR± = 



R(H ± ) - R(0) 



(8) 



where R(H^) and R(0) are the resistances for two mu- 
tually perpendicular, out-of-plane and in-plane config- 
urations of saturated magnetization, i. e., for magne- 
tization along [001] and [100] crystallographic axis, re- 
spectively. It is worth noting that under the presence 



of spin-orbit interaction, a relatively large change in re- 
sistance is expected when the direction of magnetization 
alternates from perpendicular to parallel in respect to 
the current, even if the effect of epitaxial strain, which 
makes the [100] and [001] directions non-equivalent, is 
disregarded. In the Boltzmann conductance regime, the 
effect is known as anisotropic magnetoresistance (AMR) , 
and has already been studied in (Ga,Mn)As experimen- 
tally and theoretically^ 
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FIG. 5: [color on-line] The bias dependence of the rela- 
tive change in the tunneling resistance in Zener-Esaki p- 
Gao.94Mno.oeAs/n-GaAs diode when magnetization is ro- 
tated out-of-the-plane, for various hole concentrations; n = 
10 19 cm -3 . The black squares are experimental points from 
Ref. [l3 

The calculated relative changes of the structure re- 
sistance for the magnetization vector flipping between 
perpendicular-to-the-plane and in-plane directions are 
shown in Fig. [S] as a function of bias for various hole con- 
centrations in Gao.94Mno.06As. For p — 3.5 x 10 20 cm~ 3 , 
the maximum of the computed TAMR effect, exceeding 
20%, is seen at small bias voltages. These results are com- 
pared with the experimental findings of Ref. 12, where 
the structure containing ferromagnetic Gao.94Mno.06As 
with Tc ~ 70 K was studied. According to the p-d Zener 
models such a value of Tc corresponds to p « 10 20 cm~ 3 . 

We see in Fig. that the theory describes correctly 
the experimental magnitude of TAMRj^ for small bias 
at both polarizations. We see also that the computed 
TAMRj^ tends to vanish with the increase of the reverse 
bias, whereas when the forward bias is assumed, TAMRj^ 
changes sign and tends to infinity. Such change of sign for 
the forward bias is also revealed experimentally, but the 
measured TAMRj_ appears to vanish for higher values of 
positive bias. This inconsistency can be explained by re- 
calling that the computed tunneling current stops to flow 
above the tunneling cutoff voltage, which is determined 
by a sum of the energy distance from the hole Fermi level 
E F to the top of the valence band in (Ga,Mn)As and the 
energy difference between the bottom of the conduction 
band and electron Fermi level Ep in GaAs. In the exper- 
iment, however, some current related to band-gap states 



appears to dominate near the cutoff voitage. In turn, 
series buik resistances, which are not taken into account 
in the calculations, may dominate at high reverse bias. 
Accordingly, standard AMR appears to contribute to the 
experimental value of TAMRj^ in this bias regime^ 

Figure shows also the TAMRj^ calculated for dif- 
ferent hole concentrations in the magnetic layer. Due 
to the strong p-d exchange and large spin splitting in 
the (Ga,Mn)As valence band, for the low value p = 
1 x 10 19 cm~ 3 all spin subbands above the Fermi en- 
ergy have the same spin polarization and thus TAMRj_ 
does not change the sign upon applying the positive bias. 
When, however, the Fermi level is very deep in the band, 
in the case of large p = 1 X 10 21 cm' 3 , different spin sub- 
bands contribute to the current for various voltages and 
TAMR^ as a function of bias changes the sign twice. 
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FIG. 6: [color on-line] The calculated bias dependence of the 
TMR effect in p-Gai_ !c Mn !c As/(GaAs)4/p-Gai- i EMn !C As tri- 
layer (with the GaAs barrier width d — 4 monolayers) for 
various hole concentrations p and Mn ions content x — 0.08 
(a) and for various values of x and p = 3.5 x 10 20 cm -3 (b). 



IV. TUNNELING MAGNETORESISTANCE 
A. Bias dependence 

In the previous section, we have used the model to con- 
sider a device with just one interface between magnetic 
(Ga,Mn)As and nonmagnetic GaAs. The typical tun- 
neling magnetoresistance (TMR) devices are, however, 
more complicated - they consist of a trilayer structure 
with two such interfaces, for instance, two magnetic p- 
type Ga^Mni-^As contacts separated by a nonmagnetic 
GaAs barrier. In such structures a strong TMR effect, 
i. e., a large difference in the resistance of the device for 
two configurations: parallel (ferromagnetic - FM) and 
the antiparallel (antiferromagnetic - AFM) alignments 
of magnetizations in the contacts, has been observedi 2 ^*^ 
The TMR value is usually described by the ratio, 



TMR: 



Rfm — -Rafm 



Ra 



FM 



(9) 



where Rfm and R-afm are the structure resistances for 
the FM and AFM configuration, respectively. Similarly 
to the spin polarization of the tunneling current in Zener- 
Esaki diode, TMR increases with the content of the mag- 
netic ions and decreases with the concentration of the 
holes in (Ga,Mn)As layers^ Unfortunately, in all exper- 
iments the observed TMR shows a rapid and hitherto 
unexplained decay with the increase of the applied bias. 
As shown in Fig. |SJ our calculations reproduce such de- 
cay. This suggests that the dependence of TMR on the 
applied bias results predominantly from the band struc- 
ture effects in this case. 

In Fig. Eta) one should note a strong dependence of 
TMR and its decay with applied voltage on the hole con- 
centration in the magnetic contacts. However, one can 
also see that the hole concentration does not influence 
very strongly the bias where the TMR reaches zero. This 
voltage, ca 0.3 V, corresponds to the valence band offset 
between Gao.92Mno.osAs and GaAs, which is determined 
by the spin splitting in the valence band of the former. 



The presented in Fig. EJb) TMR ratios for the magnetic 
contacts with various Mn content, i e., with different spin 
splitting, confirms this conclusion. These results suggest 
that the TMR and its decrease with the applied bias can 
be controlled by appropriate engineering of the band off- 
sets in the heterostructure, in particular, by a proper 
choice of the nonmagnetic barrier. To check this predic- 
tion theoretically, we replace in the calculations the GaAs 
by AlAs, which produces a higher by 0.55 V barrier for 
the holes. The results presented in Fig. [7| show that, in- 
deed, for tunneling through the AlAs barrier, the TMR 
magnitude decreases with the bias much slower than in 
the case of GaAs. Moreover, for a higher tunneling bar- 
rier one can expect also higher TMR ratios. 
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FIG. 7: [color on-line] The bias dependence of the TMR ratio 
for p-Gai_ a; Mn a; As/(AlAs)4/p-Gai_ 2 ,Mn a ,As with x = 0.08 
and various hole concentrations (a); for various values of x 
and p = 3.5 x 10 20 cm" 3 (b). 



A related behavior can be seen in the study of the 
dependence of the TMR effect on the width of the bar- 
rier, presented in Fig. [8] In agreement with experimental 
observations^ the calculated TMR drops rapidly when 
the barrier becomes wider. However, for AlAs which 
forms a higher barrier, the decrease of the TMR ratio 
with the number of barrier monolayers is much weaker. 
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FIG. 8: [color on-line] The bias dependence of the TMR 
ratio for various thicknesses d of the barrier layer in 
(a) p-Gao.92Mn .o8As/(GaAs) (i /p-Gao.92Mno.o8As; (b) p- 
Gao.92Mn .o8As/(AlAs) d /p-Gao.92Mn .o8As 
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B. Anisotropy of tunneling magnetoresistance 



The calculated dependence of TMR on the in-plane di- 
rection of the magnetization vector is shown in Fig. [5] It 
is seen that in the case of TMR the [100] and [110] mag- 
netization directions remains not equivalent, while TMR 
is identical for [110] and [110]. Thus, the Did symmetry 
is recovered if two interfaces are involved, in contrast to 
the case of spin current polarization in the Esaki-Zener 
diode, where Civ symmetry of a single zinc-blende in- 
terface led to the non-equivalence of the [110] and [TlO] 
directions, as discussed in the previous section. 

As shown in Fig. \Q the in-plane anisotropy of TMR 
depends crucially on the hole concentration in the mag- 
netic layer. For hole concentrations p in the range of 
10 20 cm" 3 the obtained anisotropy of TMR is below 10%, 
however, for low concentrations, p = 10 19 cm" 3 , it be- 
comes as strong as 250%. The reason for this behavior 
becomes clear when we look at Fig. \n]\ that shows the 
dependence of the tunneling current in the AFM config- 
uration on the in-plane wave vector for different direc- 
tions of magnetization. In Fig. llOf a). i. e., in the case of 
p = 10 19 cm" 3 the region of the Brillouin zone that takes 
part in the tunneling is strongly dependent on the mag- 
netization direction, in contrast to the results for higher 
hole concentrations presented in Fig. llOf b'). It should be 
stressed that such behavior has been obtained only for 
the AFM alignment. The calculated tunneling current in 
the FM configuration does not virtually depend on the 
direction of magnetization - only a very small difference 
between [110] and [010] direction has been found. 

A similar effect can be also noticed in the calculated de- 
pendence of the TMR ratio on strain, as shown in Fig. II II 
Here again we see that upon trigonal strain the tunneling 
current becomes anisotropic only for the AFM alignment 
of magnetization in the two magnetic contacts. 



FIG. 9: [color on-line] Dependence of TMR ratio on the 
direction of in-plane magnetization for various hole concen- 
trations in p-Gao.92Mno.o8As/(GaAs)4/p-Gao.92Mno.o8As tri- 
layer structure. 
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FIG. 10: [color on-line] Dependence of tunneling current on 
the direction of the in-plane wave vector for antiparallel con- 
figuration of magnetizations (AFM) along the [110] and [010] 
crystallographic directions, as indicated in the plots, and for 
hole concentrations 1 x 10 19 cm" 3 (a) and 3.5 x 10 20 cm" 3 
(b). 



C. Tunneling anisotropic magnetoresistence 

The results presented above show that the com- 
puted anisotropy of TMR results exclusively from the 
anisotropy of tunneling in the AFM configuration. Thus, 
this anisotropy cannot explain the in-plane TAMR 
effect observed in the (Ga,Mn)As/GaAs/Ga,Mn)As 
structures^*!* i. e., the difference in the resistance of 
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and (Ga,Mn)As/AlAs/(Ga,Mn)As structures. The re- 
sults are presented in Fig. In the case of the GaAs 
spacer, for hole concentrations of about 1 x 10 20 cm -3 a 
very small effect that weakly depends on the applied bias 
can be observed. Still, for small hole concentration we see 
a positive TAMR ± of the order of 60%. The difference 
between the out-of-plane and in-plane resistances of the 
same sign and of about 12% was experimentally observed 
in Ref. 14, As shown in Fig. 1131 again a higher barrier, 
i. e., the AlAs spacer, should enhance the TAMRj^ effect. 



b) 



FIG. 11: [color on-line] Dependence of tunneling current on 
the direction of in-plane magnetization for (a) parallel (FM) 
and (b) antiparallel (AFM) magnetization orientations in 
p-Gao.92Mno.o8As/(GaAs)4/p-Gao.92Mno.o8As trilayer struc- 
tures trigonally distorted along the [110] crystallographic axis; 
p = 3.5 x 10 20 cm" 3 . 



the structures in FM configuration between the x and 
y directions. As stated before, we obtain only a very 
small difference in the FM tunneling current between the 
[110] and [010] direction. Although this anisotropy can 
be increased to about 15% by assuming a low hole con- 
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FIG. 13: [color on-line] Bias dependence of the TAMR X 
ratio for TMR structure consisting of two magnetic p- 
Gao.92Mno.08As layers separated by (a) GaAs and (b) AlAs 



4-- c 1 1 nl9 —3 / t~v pt7!1 .1 barrier layer; p = 3.5 x 10 cm 

centration of p = 1 X 10 a cm 13 (compare Fig. I12|) . the J ' y 



depletion effects would not affect the anisotropy direc- 
tions determined by the symmetry of the structure. 
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FIG. 12: [color on-line] Dependence of tunneling current on 
the direction of in-plane magnetization structures with par- 
allel magnetization configuration (FM) for various hole con- 
centrations in p-Gao.92Mno.o8As/(GaAs)4/p-Gao.92Mno.o8As 
trilayer structure; p = 3.5 x 10 20 cm -3 . 



However, the calculated tunneling current for magne- 
tization vector perpendicular to the plane differs from 
the current calculated for the in-plane magnetization vec- 
tor even for the FM configuration. Using the TAMRj^ 
ratio defined in Eq. JHJ, we have calculated the per- 
pendicular TAMR for the (Ga,Mn)As/GaAs/(Ga,Mn)As 



V. SUMMARY 

We have developed the model of quantum transport 
in spatially modulated structures of hole-controlled di- 
luted ferromagnetic semiconductors, taking into account 
relevant features of the band structure within the tight- 
binding approximation. The model disregards disor- 
der and effects of carrier-carrier interactions, so that 
it is applicable to the carrier density range and length 
scales, where localization effects are unimportant. The 
computation results presented in this and our previous 
paper a 1 ^ 1 ^ demonstrate that many of experimentally im- 
portant effects, such as large magnitudes of both spin po- 
larization of the tunneling current in Zener-Esaki diodes 
and TMR ratio in trilayer structures can be understood 
within the proposed model. Furthermore, the theory de- 
scribes quantitatively a fast decay of the spin polarization 
of the current and TMR with the bias voltage without 
invoking inelastic processes. However, these processes to- 
gether with heating of Mn spin subsystem may become 
crucial in the highest bias regime. The detail studies of 
anisotropy effects reveal the presence of Civ symmetry in 
the magnitude of current spin polarization in the Zener- 
Esaki diode. This indicates a noticeable importance of 
inversion asymmetry terms specific to interfaces and zinc- 
blende structure in tunneling structures. These effects 
are not taken into account within the standard fcp-type 
approaches. According to our findings, if strain is not 
excessively large, the dominant anisotropy appears when 
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the direction of magnetization changes from parallel to 
perpendicular to the current, in a full analogy to AMR. 
Finally, we have used the model to predict theoretical 
conditions for improving the performance of the studied 
devices. Our results indicate that an attempt to increase 
the electron concentration in the n-GaAs layer of the 
Zener-Esaki tunnel junction should pay off in an increase 
of the spin polarization of the current. For the trilayers, 
the calculations suggest that reducing barrier thickness 
and increasing barrier height may result in higher values 
of TMR and its slower decay with the applied bias. 
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